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ABSTRACT 

We report results from a 2007 Suzaku observation of the Seyfert 1 AGN NGC 4593. The narrow 
Fe Ka emission line has a FWHM width ~ 4000 km s -1 , indicating emission from > 5000 R g . There 
is no evidence for a relativistically broadened Fe K line, consistent with the presence of a radiatively 
efficient outer disk which is truncated or transitions to an interior radiatively inefficient flow. 

The Suzaku observation caught the source in a low-flux state; comparison to a 2002 XMM-Newton 
observation indicates that the hard X-ray flux decreased by 3.6, while the Fe Ka line intensity and 
width a each roughly halved. Two model-dependent explanations for the changes in Fe Ka line profile 
are explored. In one, the Fe Ka line width has decreased from ^10000 to ^4000 km s _1 from 2002 
to 2007, suggesting that the thin disk truncation/transition radius has increased from 1000-2000 to 
>5000 i? g . However, there are indications from other compact accreting systems that such truncation 
radii tend to be associated only with accretion rates relative to Eddington much lower than that 
of NGC 4593. In the second model, the line profile in the XMM-Newton observation consists of 
a time-invariant narrow component plus a broad component originating from the inner part of the 
truncated disk (~300 R s ) which has responded to the drop in continuum flux. The Compton reflection 
component strength R is ~ 1.1, consistent with the measured Fe Ka line total equivalent width with 
an Fe abundance 1.7 times the solar value. The modest soft excess, modeled well by either thermal 
bremsstrahlung emission or by Comptonization of soft seed photons in an optical thin plasma, has 
fallen by a factor of ^20 from 2002 to 2007, ruling out emission from a region 5 lt-yr in size. 
Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies — galaxies: individual 
(NGC 4593) 



1. INTRODUCTION 

Accretion onto supermassive black holes in many Ac- 
tive Galactic Nuclei (AGN) is generally thought to pro- 
ceed via a radiatively efficient, optically thick, geometri- 
cally thin disk ("a-disk" , e.g., Shakura & Sunyaev 1973), 
as evidenced by their optical/UV continua (the "big blue 
bump"; Sun & Malkan 1989). Galactic black hole sys- 
tems also contain evidence for such a component; the 
thermal blackbody emission extends into the soft X-ray 
band. In Seyferts, the Fe Ka emission line at 6.4 keV 
line is a key tracer of the radiatively efficient circum- 
nuclear material. Narrow lines with Doppler-broadened 
FWHMs of a few thousand km s _1 are virtually ubiqui- 
tous in X-ray spectra of low-redshift Seyferts observed 
with XMM-Newton, Chandra-HETGS or Suzaku (e.g., 
Yaqoob & Padmanabhan 2004, Nandra 2006) and indi- 
cate material originating light-days or farther from the 
black hole, e.g., in the outer disk or molecular torus. Fe 
Ka line emission originating near the innermost stable 
orbit of the disk yields a broad (FWHMs ~ 0.1c) and red- 
shifted profile sculpted by general and special relativistic 
effects in the regime of strong gravity (Fabian et al. 1989, 
2002), although to accurately gauge the strength of the 
line, one must correctly model any continuum curvature 
which may be associated with ionized absorption along 
the line of sight (e.g., Reeves et al. 2004, Turner et al. 
2005). 

However, the flow in compact systems accreting at rel- 
atively low values relative to the Eddington limit may 



contain a radiatively inefficient component (a radiatively 
inefficient accretion flow or RIAF; Narayan et al. 1998, 
Quataert 2001); a configuration consisting of an inner 
RIAF flow which is surrounded by an a-disk beyond some 
transition radius (e.g., Esin et al. 1997) has been sug- 
gested for some low-luminosity AGN (Lu & Wang 2000 
and references therein). In the case of the Seyfert 1 AGN 
NGC 4593, modeling of the optical/UV continuum in- 
dicates blackbody emission from a truncated thin disk, 
with an inner radius of emission of 30 -Rg 1 (Lu & Wang 
2000). An observation with XMM-Newton in 2002 indi- 
cated a FWHM line width of 11000±1000 km s" 1 ; as- 
suming a black hole mass Mbh of 10 68 M Q , this width 
suggests an inner extent of no less than roughly 1000 R g 
(Brenneman et al. 2007, hereafter B07). No broad Fe Ka 
line has been confirmed with XMM-Newton (Reynolds et 
al. 2004, B07). These results were consistent with those 
obtained from a Chandra-HETGS observation in 2001, 
from which Yaqoob & Padmanabhan (2004) measured 
a FWHM line width of 2140± 1 g 8 , km s _1 . An earlier 
ASCA observation in 1994 indicated line emission from a 
radius of >30-fL, consistent with the thin disk truncation 
radius suggested from optical/UV continuum ("big blue 
bump") spectral energy distribution (SED) fitting (Lu & 
Wang 2000). Furthermore, the 2002 XMM-Newton ob- 
servation revealed, in addition to the Fe Ka core, line 
emission from ionized Fe, likely Fe XXVI (B07). Such 

1 1 R g = GM Bli /c 2 
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line features are relatively rare in Seyfert X-ray spectra, 
and could potentially originate in the collisionally ionized 
transition region between the inner radiatively inefficient 
flow and the thin disk (e.g., Lu & Wang 2000). 

Another diagnostic reflection component present in 
hard X-ray spectra is the Compton reflection hump peak- 
ing at ~20-30 keV, expected when the hard X-ray power- 
law continuum illuminates optically thick material. The 
strength of the Compton reflection hump R was found to 
be ~ 1 in a 1997 observation with BeppoSAX (Guainazzi 
et al. 1999); R = 1 corresponds to a slab covering 2n sr 
of the sky as seen from the illuminating X-ray source. 
The Suzaku X-ray observatory is the first since Bep- 
poSAX to provide continuous coverage from below 1 keV 
to at least 50 keV, allowing users to spectrally deconvolve 
the broadband continuum components (absorbing com- 
ponents, the primary power-law, and the Compton re- 
flection hump), but Suzaku 's lower >10 keV background 
yields a more accurate determination of R. 

In this paper, we report results from an observation 
of the nucleus of NGC 4593 with Suzaku in 2007, with 
the goals of constraining the Fe K emission profile and 
accurately determining the strength of the Compton re- 
flection hump in order to constrain the geometry of 
the circumnuclear accreting gas. As demonstrated be- 
low, Suzaku caught the source at an atypically low 2-10 
keV flux level, a factor of almost 4 lower than during 
the XMM-Newton observation. We observe significant 
changes in the Fe Ka profile between the 2002 XMM- 
Newton and 2007 Suzaku observations which may be re- 
lated to the decrease in continuum flux. We also report 
evolution in the strength of the soft excess, and, tenta- 
tively, the ionized Fe K emission. 

The rest of this paper is organized as follows: Sec- 
tion 2 describes the observations and data reduction. In 
Section 3, we present fits to the Fe K emission com- 
plex observed with Suzaku and compare the results to 
those for the 2002 XMM-Newton observation. In Section 
4, we present fits to the 0.3-76 keV broadband Suzaku 
time-averaged spectrum, and again compare the results 
to those obtained from XMM-Newton to investigate long- 
term spectral variability of the broadband emission com- 
ponents. The results are discussed in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 4593 (z = 0.00900, Strauss et al. 1992; a = 
12h39m39.4s, 5 = -5d20m39s, J2000.0 coordinates from 
the NASA Extragalactic Database) was observed with 
Suzaku from 2007 December 15 at 04:44 UT until 2007 
December 17 at 23:20 UT (observation ID 702040010). 
We used data from both the X-ray Imaging Spectrome- 
ter (XIS; Koyama et al. 2007) CCDs and the Hard X-ray 
Detector (HXD; Takahashi et al. 2007). For the HXD, 
we used only data collected with the PIN diodes; we did 
not consider GSO data in the paper due to the faintness 
of the source relative to the GSO background. The XIS- 
nominal pointing position was used. Further details of 
the Suzaku observatory are given in Mitsuda et al. (2007). 

2.1. XIS Reduction 

The XIS data used in this paper were version 2.1.6.16 
of the screened data provided by the Suzaku team; data 
collected within 436 s of passage through the South At- 
lantic Anomaly (SAA) were discarded, and data were 



selected to be at least 5° in elevation above the Earth 
rim (20° above the day-Earth rim). All XIS data were 
taken in normal clocking mode. 

The XIS consists of four CCDs, numbered 0-3. XIS-0, 
2, and 3 are front- illuminated (FI), the fourth (XIS-1) is 
back-illuminated (BI) and features an enhanced soft X- 
ray response. Use of XIS-2 was discontinued after 2006 
November (when it was damaged due to a likely micro- 
meteoroid impact). The CLEANSIS script was used to 
remove hot or flickering pixels. Source spectra were ex- 
tracted from a 3' radius centered on the source. The 
background was extracted using four circles of radius 
1.5', each located ~6' from the source. The XIS-FI CCDs 
were in 3x3 and 5x5 editmodes, for a net exposure time 
after screening of 118.8 ks per XIS. XIS-1 was also in 
3x3 and 5x5 editmodes, for a net exposure of 118.4 ks. 
Response matrices and ancillary response files (ARFs) 
were generated for each XIS independently using xissim- 
RMFGEN version 2007-05-14 and xissimarfgen version 
2008-04-05 (Ishisaki et al. 2007). The XIS and 3 source 
and background spectra were added using mathpha; re- 
sponse files and ARFs were added using ADDRMF and 
ADDARF, respectively. All spectra were binned to a min- 
imum of 50 counts bin -1 to allow use of the % 2 statistic 
(Gehrels 1986). 

To examine the accuracy of the XIS RMFs and deter- 
mine residual line width due e.g., to imperfect CTI cor- 
rection, we generated spectra for the emission lines gener- 
ated by the 55 Fe calibration sources illuminated two cor- 
ners of each XIS. Using the above response matrices and 
ARFs, and XSPEC v.ll.3.2ag, we fit the spectrum for 
each calibration source with a model consisting of three 
Gaussian components to model emission from Mn Kai, 
Ka2, and Kj3 (expected line centroids of 5.899, 5.888, and 
6.490 keV, respectively). The Mn Kai line energy cen- 
troids for the co-added XIS-FI, XIS-0, XIS-3, and for the 
XIS-BI were 5.908 ± 0.001, 5.906 ± 0.002, 5.910 ± 0.002, 
and 5.902 ± 0.003 keV, respectively; such discrepancies 
are consistent with the accuracy (~0.2% at the Mn Ka 
energy) of the energy calibration of the XIS. Fitting the 
co-added XIS-FI calibration source lines without the re- 
sponse file, we determined the FWHM energy resolution 
during the observation to be 168 eV. 

The 0.4-10 keV light curve, combined from all three 
operating XlSes and binned to 5760 s, is plotted in the 
top panel of Fig. 1; the variability amplitude F var (see 
Vaughan et al. 2003 for definition) is 16.8 ± 0.8%. 

2.2. HXD-PIN Reduction 

The HXD-PIN is a non-imaging instrument with a 34' 
square FWHM field of view. 12-76 keV PIN source spec- 
tra were extracted from version 2.1.6.16 event files pro- 
vided by the HXD instrument team. PIN background 
count rates are variable and strongly depend on the 
time since SAA passage (Kokubun et al. 2007), data 
were selected according to the following criteria: at least 
500 s since SAA passage, geomagnetic cutoff rigidity 
(COR) > 8 GV, and day- and night-Earth elevation an- 
gles each >5°. Instrumental (non- X-ray) background 
spectra for the PIN were provided by the HXD Team, 
who generated the background using the calibrated GSO 
data for the background monitor ("tuned" background, 
with METHOD=LCFITDT). The current accuracy of 
the PIN non- X-ray background (NXB) model for a 1 day 
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observation is <1.5% (la peak-to-peak residuals) below 
^50 keV (Fukazawa et al. 2009). Both the source and 
background spectra were generated with identical good 
time intervals, and the source exposure was corrected 
for instrument dead time (a 6.8% effect; the background 
event files were already dcadtimc-corrected by the HXD 
Team). This yielded a good exposure time of 90.2 ks. 
To model the contribution to the total background from 
the Cosmic X-ray Background (CXB), a spectrum of 
the form 9.0 x 10" 9 (£;/3kcV)" - 29 exp(-£'/40keV) erg 
cm~ 2 s _1 sr _1 keV -1 (Gruber et al. 1999) was used (see 
the Suzaku ABC Guide for further details). The total 
(X-ray plus particle) background 12-76 keV flux was 8.6 
x 10~ 10 erg cm~ 2 s _1 . 

The source spectrum was binned to a minimum 
of 400 counts bin -1 . We used the response file 
ae_hxd_pinxinome4_20080129.rsp. The mean 12-76 keV 
net source flux and count rate were 4.9 x 10~ n erg cm~ 2 
s _1 and 0.077 ct s _1 , respectively. The 12-76 keV net 
source light curve, binned to 11520 s (two satellite or- 
bits), is plotted in the bottom panel of Fig. 1; error bars 
include a ~2% systematic uncertainty (Fukazawa et al. 
2009). The variability amplitude F var was < 13.8%. 

There is a possible contaminating X-ray source, 1WGA 
J1239.7-0526, located about 5' south of' NGC 4593, and 
thus in the field of view of the PIN. However, as demon- 
strated in the Appendix, we estimate that in the 12-76 
keV band, this source contributes only 7x 10~ 13 erg cm~ 2 
s- 1 to the observed 12-76 keV spectrum of NGC 4593 
(about 1.5% of the net source flux), and can be safely 
ignored. 

3. MODEL FITS TO THE FE K BANDPASS SPECTRUM 

All spectral fitting in this paper used XSPEC 
v.ll.3.2ag. The abundances of Lodders (2003) were used. 
In all fits below, a neutral Galactic column of 1.89 x 
10 20 cm~ 2 was included (Kalberla ct al. 2005). Unless 
otherwise stated, all errors below correspond to A\ 2 — 
2.71 (90% confidence level for one interesting parameter 
when the errors are symmetric) with the XIS-BI/XIS- 
FI normalization left free (best-fit values were usually 
close to unity) and the PLN/XIS normalization set at 
1.16 (Maeda et al. 2008). 

Previous soft X-ray observations of NGC 4593 have 
indicated absorption features attributed to an outflow- 
ing, warm absorber system. Using data from a 108 
ks Chandra-LETGS observation in 2001 February and 
a 10 ks XMM-Newton observation in 2000 July, Steen- 
brugge et al. (2003) modeled two zones of ionized ab- 
sorption along the line of sight: a highly-ionized zone 
with log£ = 2.61±0.09 erg cm s _1 and a column density 
Ah = 1.6 ± 0.4 x 10 21 cm~ 2 , plus a more lowly-ionized 
zone with ionization parameter log£ near 0.5 erg cm s _1 
and Ah near 6 x 10 19 cm~ 2 ; £ = Li on n~ 1 r~ 2 , where L; on 
is defined as the isotropic 1-1000 Ryd ionizing contin- 
uum luminosity, n c is the electron number density, and 
r is the distance from the central continuum source to 
the absorbing gas. Both zones were measured to be out- 
flowing at -400 km s _1 relative to systemic. From an 
80 ks observation with Chandra-HETGS in 2001 June, 
McKernan et al. (2003) confirmed the highly-ionized ab- 
sorber, deriving A H = 5.37±Q;yg x 10 21 cm~ 2 and log£ 
= 2.521q;q4 erg cm s _1 , although they derived a lower 



outflow velocity relative to systemic. The high-column, 
high-ionization component manifests mainly itself via H- 
and H-likc lines, the most prominent ones being N VII, 
O VIII, Ne IX, Ne X, Mg XII, Si XIII, and Si XIV, as 
well as an absorption edge due to O VIII, There is also 
line absorption due to Fe L XX-XXV which, in CCD 
resolution spectra, yields an absorption trough near 0.9 
keV and a series of blended edges which yield a trough 
from ~1.3 to 2.2 keV. The low-column, low-ionization 
component manifests itself via O V and O VI absorption 
and is not expected to contribute strongly to CCD res- 
olution spectra above 0.2 keV. McKernan et al. (2003) 
also reported a neutral Fe L3 edge at 707 eV due to dust 
along the line of sight in the host galaxy of NGC 4593 
(this feature was also detected by B07 in their analysis of 
the XMM-Newton pn spectrum). A moderately strong, 
steep soft excess has been known since an EXOSAT ob- 
servation by Pounds & Turner (1988) and has been phe- 
nomenologically modeled using thermal bremsstrahlung 
and Comptonization models (e.g., B07). 

To fit the Fe K bandpass, we restricted ourselves to 
the >4-11.5keV band of the co-added XIS-FI spectrum, 
where warm and neutral absorption and the soft excess 
are not expected to have an impact. Data/model resid- 
uals to a simple power-law are plotted in Figure 2a, and 
reveal that in addition to the clear Fe Ka emission line 
at 6.4 keV, there is emission near 7.0 keV (all photon en- 
ergies are rest frame unless otherwise stated). A model 
consisting of a power-law plus a Gaussian component to 
model the Fe Ka line yields x 2 /dof —360. 8/318; residuals 
are plotted in Fig. 2b. We added a Gaussian with energy 
centroid fixed at 7.056 keV to model Fe K/3 emission, 
with the normalization held at 0.13 times that of the Fe 
Ka emission line (i.e., assuming K/3 emission from neu- 
tral Fe; untying this parameter did not yield a significant 
improvement in this or any subsequent fit). \ 2 /dof fell 
to 332.3/318. Data/model residuals are plotted in Fig. 
2c. To model the effects of an Fe K edge at 7.11 keV 
due to Compton reflection, we added a pexrav compo- 
nent (Magdziarz & Zdziarski 1995), with an inclination 
angle of 30° assumed, and a power-law cutoff set at 500 
keV. The reflection strength R was kept fixed at 1.08 and 
the Fe abundance was kept at solar, as per the best-fit 
model to the 0.3-76 keV Suzaku spectrum (see Section 
4). x 2 /dof fell to 320.6/318. This is our best-fit model 
to the Fe K bandpass spectrum; best-fit parameters are 
given in Table 1. 

Data/model residuals are plotted in Fig. 2d. There 
appear to be some moderate (~7%) positive residuals at 
6.7 keV, the energy of Fe XXV. However, there was no 
significant improvement to the fit when we added another 
Gaussian component with energy centroid fixed at 6.70 
keV and width a fixed at 1 eV. \ 2 fell by only 2.5 for 
one less dof, significant at only the 88% confidence level 
according to an F-test (and at 87.5% according to Monte 
Carlo simulations.) 2 The upper limit on the intensity 

2 As pointed out by Protassov et al. (2002), it is improper to 
use the F-test to determine the statistical significance associated 
with comparing a model which contains an emission or absorption 
line to a "null hypothesis" model where the line does not exist. 
This is because the F-test cannot be used when the null values of 
one parameter lie on the boundary of possible parameter values. In 
the case of any emission or absorption line, the boundary condition 
of zero line flux occurs for the null hypothesis model. When the 
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(equivalent width, EW) to Fe XXV emission was 2.8 x 
10~ 6 ph cm- 2 s" 1 (19 eV). 

B07 reported detection of an emission line near 6.9 
keV, likely associated with Fe XXVI. We added to our 
model a narrow (width cr=l eV) Gaussian component 
with energy centroid fixed at 6.96 keV, but no there was 
no significant improvement to the fit, as |A% 2 | was < 
2.71 for one less dof. The upper limit on the intensity 
(EW) to Fe XXVI emission was 1.8 x 1CT 6 ph cm" 2 s"' 1 
(17 eV). 

The Fe Ka line's width a is 41±} 2 eV (41±} 2 eV after 
subtracting in quadrature the instrumental broadening 
of <9 eV based on the fits to the 55 Fe calibration sources 
in Section 2.1). The line EW is 255±19 eV. The best-fit 
energy centroid of the Fe Ka line in the co-added spec- 
trum is 6.421±0.007 keV, or 6.412±0.007 keV after ac- 
counting for the small gain offset as indicated by the 55 Fe 
calibration line. The best-fit line energies in separate fits 
to the individual XIS-0 and XIS-3, and XIS-1 spectra are 
6.422±0.009, 6.408±0.009 keV, and 6.401±0.014keV, re- 
spectively, after accounting for the gain offsets in each of 
those detectors. The small implied offset from 6.40 keV 
in the co-added spectrum is likely not real, given that 
XIS-3 and XIS-1 each indicate emission consistent with 
neutral Fe, and there may be a slight energy scale prob- 
lem associated with XIS under investigation as of this 
writing. We also note that the Fe Ka line width in each 
XIS is consistent with that from the co-added spectrum. 

We tested for the presence of a Compton shoulder, ex- 
pected if the bulk of the narrow Fe Ka line (hereafter 
also referred to as the "core") originates in Compton- 
thick material. We modeled such a component using a 
Gaussian component with energy centroid fixed at 6.32 
keV and width a left free. 3 We found an upper limit 
to the intensity (EW) to Compton shoulder emission of 
7 x 1(T 6 ph cm~ 2 s- 1 (44 eV), or about 24% that of the 
Fe Ka core. 

F-test is used in this manner, the resulting statistical significance 
is at best only an "estimate" of the true stastical significance. In 
this paper, we include the results from the F-test, but we also 
include results using Monte Carlo simulations, described in detail 
in Section 3.3 of Porquct ct al. (2004). In these simulations, one 
searches over the possible energy ranges where a feature may be 
expected to be detected, fitting an unresolved Gaussian at multiple 
energy bins to test if fitting spurious features due to photon noise 
can yield a reduction in \ 2 as large as the reduction associated 
with including the emission/absorption line in the model. In the 
case of a line with an expected observed energy, e.g., 6.70 keV in 
the case of an expected Fe XXV emission line, we searched over an 
energy range which was centered on the line expected energy and 
had a width ± the FWHM resolution of the instrument. Through- 
out this paper, we list the statistical significances determined from 
Monte Carlo simulations in parentheses following the significances 
estimated using the F-test. 

3 The sum of first scatterings of Fe Ka photons forms a "shelf" - 
like feature extending redward from 6.400 keV, with the maximum 
downward energy shift due to a single Compton scattering being 
0.156 keV; the shape of the "shelf" and the intensity relative to 
that of the Fe Ka core arc dependent upon geometry, column den- 
sity of the reflecting material, and viewing angle (see e.g., Murphy 
& Yaqoob 2009; Watanabe et al. 2003; Matt 2002). Given the 
quality of the current data set and the energy resolution of the 
XIS, a Gaussian emission component can be an adequate substi- 
tute for a first-scattering Compton shoulder if the energy centroid 
is constrained to >6.244 keV (ideally, centered near 6.32 keV) and 
the width a is ~0.1-0.2 keV. 



3.1. Comparison with the 2002 XMM-Newton 
Observation 

To investigate possible evolution of the Fe K emission 
complex, we compared the results obtained with Suzaku 
to those obtained in a re-analysis of the 76 ks observation 
in 2002 June by the XMM-Newton (ObsID 0059830101). 
We downloaded Observation Data Files (ODF) events 
data and reprocessed events for the European Imaging 
Photon Camera (EPIC) pn detector using XMM Science 
Analysis Software v. 7. 1.0. Further details of the obser- 
vation are given in B07. Source and background spectra 
were extracting using standard analysis methods for a 
good exposure time of 53 ks (the pn had been in small 
window mode to prevent pile- up). The source spectrum 
was grouped to a minimum of 40 ct bin -1 . 

We fit the 4-11.5 keV pn spectrum with a model con- 
sisting of a power-law component plus Gaussian com- 
ponents to model emission from Fe I Ka and K(3, and 
Fe XXVI. We also included a Compton reflection com- 
ponent using pexrav with R fixed at 1.08. 4 It was not 
significant to include an edge at 7.11 keV (optical depth r 
< 0.03), so the choice of R is somewhat arbitrary given 
that there was no simultaneous > 10 keV energy cov- 
erage. However, the depth of the Fe K edge does not 
significantly affect the conclusions, and so we include 
this component for completeness. The best-fit param- 
eters are listed in Table 1 for comparison with the 2007 
Suzaku results; parameters are consistent with those ob- 
tained by B07 for the Fe Ka and the Fe XXVI lines. Fig. 
3a shows data/model residuals to a model consisting of 
only a power-law and a Compton reflection hump for 
both the XMM-Newton and Suzaku Fe K bandpass spec- 
tra; Fig. 4 shows contour plots of line intensity versus 
centroid energy to this model indicating the significance 
of including an unresolved Gaussian component. Fig. 3b 
shows the data/model residuals when the Fe Ka and Fe 
Kfi lines are additionally modeled. 

We also tested for a possible Compton shoulder, mod- 
eled the same way as in Section 3.1; the upper limit to 
Compton shoulder emission intensity (EW) was 1.1 x 
10~ 5 ph cm" 2 s- 1 (22 eV), or 23% that of the Fe Ka 
core. 

Comparing the results from the 2007 Suzaku and 2002 
XMM-Newton observations, several things are apparent. 
While the 4-10 keV flux has decreased from 2.47 x 10~ n 
erg cm -2 s" 1 in 2002 to 0.69 x 10~ n erg cm" 2 s" 1 in 
2007 (a factor of 3.6), the total intensity of the Fe Ka line 
has decreased from 4.79 ± 0.48 x 10~ 5 (2002) to 2.93 ± 
0.22 x 10~ 5 ph cm" 2 s" 1 (2007), a factor of only 1.6. 
The Fe Ka line EW thus has increased from 114 ±11 eV 

4 Fixing R at 1.08 follows the assumption that the Compton 
reflection component has tracked the continuum flux over 5 years. 
If the absolute normalization of the Compton hump tracked the 
intensity of the Fe Ka line from 2002 to 2007, then we would expect 
R to increase by a factor of 2.24—2.42 (consistent with tracking the 
EW of the Fe line). In the 2002/ XMM-Newton spectrum, if R were 
really 2.42, then we would observe a significant Fc K edge at 7.1 
keV. With the lack of > 10 keV coverage, the value of R is poorly 
constrained, but fixing R at 2.42 (or any value > 1), with solar Fc 
abundance, yielded edge-like residuals in emission above 7.1 keV. 
Thawing the Fe abundance did not yield data/model residuals as 
good as when R had been fixed at 1.08. The exact choice of R 
in the fits is thus somewhat arbitrary, but does not significantly 
affect our conclusions in the paper regarding the evolution in the 
observed width of the Fc Ka line or in the soft excess. 
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(2002) to 255±19 eV (2007). In this model where the Fe 
Ka line is modeled with a single Gaussian component, 
the width has decreased from 87±17 eV (2002) to 4l±[ 2 6 
eV (2007). The two line widths arc inconsistent up to 
the Ax 2 = 9.5 (>99.7%, or >3cr) confidence level. Given 
that the XIS and the pn have similar energy resolution 
near 6 keV, the evolution in line width is likely intrinsic 
to the source. 

Finally, there is highly tentative evidence that the ion- 
ized emission blucward of 6.4 keV has evolved, too. The 
lower limit to the intensity of the Fe XXVI line in the pn 
spectrum and the upper limit to Fe XXVI emission (as- 
suming an unresolved line in each case) are just barely 
consistent at the 90% confidence level Figs. 3b and 4 hint 
at a marginal change in line intensity as well. We cau- 
tion, however, that, given the energy resolution and the 
quality of the spectrum, there is likely blending between 
the Fe XXVI and Fe Kj3 lines and Fe K edge at 7.11 keV, 
and the width of the Fe XXVI line is not constrained. 

3.2. Dual-Gaussian Fit to the Fe Ka Core 

As an alternative, we considered the possibility that 
the Fe Ka core in NGC 4593 consists of two Gaussian 
components, one narrow (width a ~ 40 eV) and one rel- 
atively broader, with a ~ 100 eV. The narrow component 
is assumed to be present in both spectra. The broad com- 
ponent is assumed to respond to the drop in continuum 
flux and assumed to be detected only in the 2002 XMM- 
Newton spectrum and to be too faint to be detected in 
the 2007 Suzaku spectrum. 

We first fit the XMM-Newton spectrum with this 
"dual-Gaussian" model. The narrow Gaussian compo- 
nent's width (Ti was fixed at 41 eV, the best-fit value 
from the Suzaku fit, and the broader Gaussian compo- 
nent's width (72 was left free. The intensity of the narrow 
Gaussian, I\ , was initially kept fixed at the best-fit value 
from the Suzaku spectrum, 2.93 x 10 -5 ph cm -2 s _1 , as 
might be expected if this component has not responded 
to the drop in continuum flux in the Suzaku spectrum; 
thawing this parameter did not yield a significant im- 
provement in the fit and constraints on 7i were poor. 
The intensity of the broad line, I 2 , was left as a free pa- 
rameter. Both Gaussian energy centroids were fixed at 
6.40 keV. In the best-fit model, x 2 /dof = 1106.3/1092, 

0-2 = 177ig eV, h = 2.18±g;|2 x 10~ 5 ph cm" 2 s" 1 . 
The EWs of the narrow and broad components in the 
model were 64 eV and 47^ 7 eV, respectively, i.e., both 
components comprise roughly equal fractions of the total 
line flux. 

Returning to the Suzaku spectrum, we added to our 
best-fit model a Gaussian component with width fixed 
at 177 eV. This did not yield a significant improvement 
in the fit; the upper limit to the intensity (EW) was 
1.0 x 10~ 5 ph cm" 2 s" 1 (65 eV). 

4. THE BROADBAND SUZAKU SPECTRUM 

We fit the XIS-FI data from 0.4-11.5 keV, the XIS-BI 
data from 0.3-10.5 keV, and the PIN data from 12-76 
keV. Due to calibration uncertainties associated with the 
instrumental Si K edge, data from 1.75-1.83 keV were 
ignored in the XIS-FI and from 1.75-1.87 keV in the 
XIS-BI data. All the data are plotted in Figure 5a. 

Data/model residuals to a simple power-law compo- 



nent absorbed only by the Galactic column are shown in 
Figure 5b and indicate spectral features reported in pre- 
vious studies: an absorption feature near 0.9 keV, the 
moderate soft excess at lower energies, and an excess > 
12 keV which is likely the Compton reflection hump. 

Our "initial guess" model was based on the best fit 
to the XMM-Newton spectrum by B07. It consisted 
of the following components: a power-law; a Comp- 
ton reflection hump modeled with PEXRAV, with abun- 
dances set to solar values, an inclination angle of 30° 
assumed, and a power-law cutoff set at 500 keV; and 
a soft excess phenomenologically parametrized as ther- 
mal bremsstrahlung emission, using the model zbrem in 
xspec. We modeled one zone of ionized absorption with 
an xstar v.2.11 table which assumed solar abundances, 
a turbulent velocity of 200 km s -1 and an input photon 
spectrum of 2. The absorber's outflow velocity relative to 
systemic was kept fixed at -400 km s -1 (Strccnbrugge et 
al. 2003) The best-fit model had x 2 /dof = 1285.4/1110; 
residuals are plotted in Figure 5c. 

We then added components to model features reported 
in previous observations. We added narrow (width a = 
1 eV) Gaussian components to model line emission from 
O VII / (energy centroid fixed at 561 eV) and Ne IX / 
(energy centroid fixed at 905 eV), features reported us- 
ing Cfcarcdra-HETGS by McKernan et al. (2003); x 2 /dof 
dropped to 1279.8/1109 and then 1272.6/1108, respec- 
tively. We next added an Fe L 3 edge with energy fixed 
at 707 eV, reported by both McKernan et al. (2003) and 
B07; x 2 /dof fell to 1262.1/1107. A final improvement in 
the fit was achieved by adding a third narrow Gaussian 
near 1.33 keV, the energy of He-like Mg, where there still 
remained ^+3 to +4% residuals. We assumed emission 
from Mg XI / and fixed the Gaussian energy centroid 
at 1.331 keV; x 2 /dof fell to 1254.5/1106. Data/model 
residuals are plotted in Figure 5d. Best-fit model pa- 
rameters are listed in Tabic 2. An additional layer of 
absorption, either neutral or ionized, was not required in 
our model. It was not significant to thaw the Fe abun- 
dance of the PEXRAV component from solar values. 

In the best-fit model, the power-law component had 
T = 1.656^0 025 j an d the warm absorber was modeled 
using a column density Nh,wa = 2.8 ± 0.8 x 10 21 cm -2 
and an ionization parameter of log£=2.39 ± 0.17 erg cm 
s" 1 . The best-fit value of R was 1.08±0.20, and Fig. 6 
shows a contour plot of R as a function of T. The uncer- 
tainties here are statistical only, and calculated assuming 
the PIN/XIS-FI normalization is fixed at 1.16. 5 Assum- 
ing ±1.5% systematics in the absolute background flux, 
the systematic uncertainty on R is an additional ±0.15 
(not plotted in Fig. 6). 

We explored alternate parametrizations of the soft ex- 
cess. Modeling the soft excess as a simple power-law 
yielded a good fit, with x 2 /dof = 1257.3/1106 for best- 
fit values of T = 3.6 ±1.3 and power-law normalization 

at 1 keV of 2.3ti 4 9 7 x 10 ~ 5 P h cm ~ 2 s_1 keV" 1 . A good 
fit was also achieved assuming thermal Comptonization 
emission, modeled using COMPTT in xspec (Titarchuk 
1994). A spherical geometry was assumed, as was an 

5 When the PIN/XIS-FI normalization was allowed to vary, it 
went to a value above 1.5, which is highly inconsistent with Suzaku 
calibration. The PIN/XIS-FI normalization was thus kept fixed at 
1.16 in all fits. 
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input soft photon temperature of 50 eV. In the best-fit 
model, x 2 /dof was 1255.8/1105 for a plasma tempera- 
ture of 53^3^ keV and optical depth <0.28. In the best- 
fit model for each cases, data/model residuals and the 
0.4-1.0 keV flux of the soft excess component were vir- 
tually identical to those obtained for the best-fit model 
assuming thermal brcmsstrahlung emission. All other 
model parameters were consistent with those obtained 
when thermal bremsstrahlung emission was modeled. 

We next tested if emission from a relativistically broad- 
ened Fe K line was required. We used a LAOR compo- 
nent in XSPEC (Laor 1991). The inner radius, outer ra- 
dius, disk inclination, and emissivity index were initially 
kept fixed at 1.5 R g , 400 R g , 30° and -3, respectively. 
Similar to Reynolds et al. (2004) for the 2002 XMM- 
Newton spectrum, no significant improvement to the fit 
was found even when thawing the inner radius or disk 
inclination {\ 2 dropped by less than 2). The upper limit 
to the broad Fe K line intensity (EW) was 4 x 10~ 6 ph 
cm" 2 s- 1 (40 cV). 

We next explored if the soft excess and Compton re- 
flection component could be modeled in a self-consistent 
manner by assuming reflection from an ionized surface. 
We used the table model reflion.mod (Ross & Fabian 
2005), and set the input photon index equal to that for 
the primary power-law. Initially, the Fe abundance rela- 
tive to solar values Zp c was kept fixed at 1.0. In the best- 
fit model, x 2 /dof was 1368/1106 for £ near 1500 erg cm 
s _1 , with poor data/model residuals. The fit was able to 
model correctly most of the soft excess flux, but under- 
predicted the strength of the Compton reflection hump, 
even as the best-fit value of T went to ~1.0. Allowing 
Z-p e to vary did not improve the fit. We next modeled 
relativistic blurring of the ionized reflected emission with 
kdblur. We kept the outer radius of blurring fixed at 
400 R g and allowed the inner radius to vary from 1.24 to 
200 R g . However, no further significant improvement to 
the fit was found for any value of the inner radius tested. 
Ionized reflection models are not discussed further. 

4.1. A Re-analysis of the 2002 XMM-Newton pn 
Broadband Spectrum 

In the best-fit model, the observed (absorbed) 2-10 
keV flux was 1.03 x 10~ n erg cm -2 s -1 , a factor of 3.8 
lower than that reported by B07 for the 2002 XMM- 
Newton observation. This is also a factor of 3.8 lower 
than the long-term average 2-10 keV flux as deter- 
mined by continuous Rossi X-ray Timing Explorer Pro- 
portional Counter Array (RXTE-PCA) monitoring from 
2004 February to 2007 December 6 (Figure 7; see also 
Summons ct al., in prep.). 

To explore long-term variability of the spectral com- 
ponents in a model-dependent fashion, we tried to apply 
our best-fit Suzaku model to the 0.2-12 keV 2002 XMM- 
Newton pn spectrum, while allowing as few parameters 
as possible to change between models to both spectra. 
As per Section 3.2, we allowed the Fe Ka line energy 

6 We downloaded all the public archive data for NGC 4593 and 
created a 2—10 keV light curve, binned for each observation, stan- 
dard extraction methods; the reader is referred to e.g., Markowitz 
ct al. (2003) for details on light curve extraction and to Edelson & 
Nandra (1999) for details regarding on PCA background subtrac- 
tion. 



centroid, width, and intensity to vary, and we added a 
narrow Gaussian component to model Fe XXVI line emis- 
sion. We modeled the warm absorber using the same 
xstar table model as for the Suzaku spectrum. 

We found that, in order to achieve a good fit and 
good data/model residuals (< ±2 — 3%, it was nec- 
essary and sufficient to thaw the power-law normaliza- 
tion, power-law photon-index, bremsstrahlung normal- 
ization and temperature, Fe L3 edge depth t FeL3 > and 
warm absorber column density A/h,wa- Best-fit values 
for these parameters are listed in Table 3, alongside the 
corresponding values from the Suzaku fit. All other pa- 
rameters were kept frozen at the values listed in Ta- 
ble 2. The best-fit model had x 2 / rf o/=1793. 2/1756; 
data/model residuals are plotted in Fig. 8. We found 
that it was not necessary to include in the model a sec- 
ond, lower-ionization warm absorber along the line of 
sight. 

As reported by B07, the total absorbed 0.5-10 keV flux 
during the 2002 XMM-Newton observation was 6.74 x 
10~ n erg cm -2 s _1 , a factor of 4.7 higher than the to- 
tal absorbed 0.5-10 keV flux during the 2007 Suzaku 
observation. However, the broadband spectral changes 
between the two spectra cannot be explained solely by 
changes in power-law normalization and photon index; 
the best-fit model to the pn spectrum in this case had 
X 2 /dof ~ 4 and unacceptable data/model residuals, par- 
ticularly below 2 keV. Our model fits indicate that the 
difference in Nu,wa between the two observations is not 
highly statistically significant, as the column densities 
are consistent at the A% 2 ~ 6 confidence level. (The 
variation in TF e L3 is not likely real, either; there is some 
mild degeneracy between Tp C L3 and -/Vr\wa, and the val- 
ues of t FeL3 are consistent at the ~99% confidence level 
(Ax 2 = 7). If this component is associated with dust in 
the host galaxy along the line of sight at least several pc 
from the black hole, then we would not expect any tem- 
poral variation in tf c l3-) Evolution in the normalization 
of the soft excess must contribute to the change in soft 
X-ray flux between the two observations. For comparison 
purposes, we now focus on the 0.4-2.0 keV unabsorbed 
flux, -Fb.4_2,unabs) of this component. 

The normalization of the soft excess in the Suzaku 
spectrum is not well constrained due to its faintness. 
Assuming thermal bremsstrahlung emission, we find 
-Fo.4-2,unabs = SJt^ 2 x 10~ 14 erg cm" 2 s" 1 for this com- 
ponent in the Suzaku spectrum (uncertainties are based 
on the error on the normalization of this component). 
The value of i 7 b.4-2,unabs of this component in the pn 
spectrum was 4.72 ±0.19 x 10~ 12 erg cm~ 2 s _1 . We con- 
clude that, in the context of modeling the soft excess as 
thermal bremsstrahlung emission, the normalization has 
dropped by at least a factor of 20. Unfolded model spec- 
tra for our best-fit models to the pn and Suzaku spectra 
assuming thermal bremsstrahlung emission are shown in 
Figs. 9 and 10, respectively. 

In the XMM-Newton spectrum, we substituted the 
thermal bremsstrahlung component with a thermal 
Comptonization component, again using COMPTT, 
and keeping the geometry parameter, input soft photon 
temperature, and plasma temperature fixed at the best- 
fit values found for the Suzaku spectrum; only the opti- 
cal depth and normalization were allowed to vary. This 
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yielded a similar good fit, with % 2 /dof = 1922.6/1755; 
best-fit values for the free parameters are listed in Ta- 
ble 3. Fo.4-2. unabs of the Comptonization component 
was 1.54+o'.i3 x 10~ n erg cm' 2 s _1 (flux uncertainty is 
from the uncertainty on the normalization of this com- 
ponent). This flux is a factor of at least 20 higher than 
the corresponding flux from the best-fit Suzaku model, 
8.2+J 4 7 6 x 1(T 14 erg cm" 2 s" 1 . 

5. DISCUSSION AND CONCLUSIONS 
5.1. Summary of Observational Results 

We have presented results from a Suzaku observation 
of the nucleus of the Scyfcrt AGN NGC 4593 in 2007 
December, and we compare our spectral fits for both the 
Fe K bandpass and the broadband X-ray spectrum with 
those obtained from a 2002 XMM-Newton EPIC-pn ob- 
servation. Suzaku caught the source at a relatively low 
X-ray flux level: the 2-10 keV continuum flux during the 
Suzaku observation was a factor of 3.8 lower. 

The Fe Ka line intensity has dropped by a factor of 
1.7, suggesting that roughly half of the total line flux 
has responded to the drop in continuum flux. One of 
our main results is that the Fe Ka line is significantly 
more narrow in the Suzaku observation. Modeling the 
line as a single Gaussian, we find that the width a has 
dropped from 87±17 eV in 2002 to 41± 1 2 eV in 2007. We 
also modeled the line using a dual-Gaussian model com- 
posed of relatively narrow and broad lines. The former 
dominates the Suzaku profile and is assumed to be time- 
invariant; in the XMM-Newton spectrum, both lines are 
modeled to have roughly equal intensity and the broad 
component has a width a = 177+52 eV- There is highly 
tentative evidence for the Fe XXVI emission line at 6.96 
keV to have dropped in intensity from 2002 to 2007, as- 
suming an intrinsically narrow (unresolved) line: in the 
Suzaku observation, the line is not significantly detected 
(EW < 17 eV). 

In our broadband fits to the 0.3-76 keV spectrum, the 
primary power-law component, commonly attributed to 
inverse Comptonization of soft seed photons in a hot 
corona (e.g., Haardt et al. 1994), was relatively flat, with 
r near 1.65. The Compton reflection component had a 
relative strength R ~ 1.08. We also model a modest 
soft excess using both thermal bremsstrahlung emission 
and thermal Comptonization of soft seed photons, sim- 
ilar to B07, and we obtain similar results. Importantly, 
we find the soft excess has dropped in flux by a factor of 
at least ^20 between the XMM-Newton and Suzaku ob- 
servations. We model one zone of absorption along the 
line of sight, the previously seen highly- ionized (log£ ~ 
2.5) zone, with a column density similar to that obtained 
by McKernan et al. (2003) and Steenbrugge et al. (2003). 
There is no strong evidence for evolution of the warm ab- 
sorber between the two observations. A relativistically 
broadened Fe Ka line was not detected in the Suzaku 
spectrum; Reynolds et al. (2004) demonstrated a similar 
result in the XMM-Newton spectrum. 

5.2. Tracing the Truncated Disk with the Fe Ka Line 

We explore two (model-dependent) scenarios to cor- 
relate the changes in Fe line intensity and profile with 
the observed drop in continuum flux. In the model 
where a single Gaussian was used to describe the Fe 



Ka profile, the width a in the 2007 Suzaku observation 
was 41^16 eV, which corresponds to a FWHM velocity 

vfwhm of 4420+^ygQ km s _1 . This is roughly commensu- 
rate with the optical broad emission line width: Peterson 
et al. (2004) reported FWHM Ha and H/3 line widths of 
3399 ± 196 and 3769 ± 862 km s" 1 , respectively. 7 As- 
suming that the line originates in gas that is in virialized 
orbit around the black hole, we can estimate the distance 
r from the black hole to the line-emitting gas. Assum- 
ing that the velocity dispersion is related to Ufwhm as 
<v 2 > = §w|wHM (Nctzcr 1990), we use GM BH = rv 2 . 
We use a black hole mass Mbh of 6.6 x 10 6 M Q , an es- 
timate based on the relation between Mbh and stellar 
velocity dispersion in Seyferts (Nelson et al. 2004). The 
best-fit reverberation mapping estimate from Peterson et 
al. (2004), 5.4 x 10 6 M Q , is consistent with this estimate. 

We find r = 6.0±2° 4 2 x 10 13 m, or 2.3±?;g lt-days. As 1 
i? g corresponds to 1 x 10 10 m for the black hole mass used, 

r = 6000^2500° R g - We cannot of course rule out contri- 
bution from an even more narrow Gaussian component 
originating in even more distant material. In the 2002 
XMM-Newton observation, the corresponding measured 
line width (we use our best-fit value of a = 87±17 eV) 
corresponds to a value of r = l-3±g;J x 10 13 m = 0.50+^ 

lt-days, or l350±H° R g (see also B07). B07 also used the 
lack of observed variability in the Fe line flux during the 
XMM-Newton observation to constrain the light-crossing 
time for the line-emitting gas to be at least 2000 i? g . 

One possible explanation to explain the change in Fe 
Ka line profile, insofar as it traces the geometrically thin, 
radiatively efficient disk, is that the innermost radius of 
the thin disk has increased over 5 years. A common 
model for accretion flows incorporating truncated thin 
disks is one where the thin disk transitions to an inner 
RIAF as the flow crosses a certain transition radius r t 
(Esin et al. 1997); a commonly invoked type of RIAF is 
an advected-dominated accretion flow (ADAF), wherein 
the disk is optically thin and geometrically thick (e.g., 
Narayan & Yi 1995). The largest width observed for 
the Fe Ka line thus could indicate r t . In this model, 
r t is expected to increase, and more of the inner thin 
disk evaporates, as the accretion rate relative to Edding- 
ton, m, decreases in a given object. Supporting evidence 
for this comes from timing observations of black hole 
X-ray Binary systems during outburst decay: character- 
istic temporal frequencies in the power spectral density 
function (PSD), such as peaks of Lorentzian components 
and/or quasi-periodic oscillations, migrate towards lower 
temporal frequencies as m decreases and the source lumi- 
nosity fades, as the source evolves through the low/hard 
spectral state into quiescence (e.g., Axelsson et al. 2005, 
Belloni et al. 2005). In addition, the temperature and 
flux of the soft, thermally emitted component have been 
seen to decrease with rh in many sources (e.g., Gierlihski, 
Done & Page 2008). 

However, the predicted relationship between m and r t 
remains unclear. Yuan & Narayan (2004) empirically 

7 Peterson et al. (2004) reported that the continuum-line lag 
results were poor. The Ha lag was reported as 3.2^'^ lt-days 
but Peterson et al. (2004) recommended caution. The H/3 lag was 
reportedly "completely unreliable." 
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derive that compact sources accreting at m near 10~ 2 , 
10~~ 4 , and 10~~( 6 ~ 7 ) may be associated with values of 
r t near 10~( 1_2 \ 10~( 2 ~ 3 \ and 10~( 4 ~ 5 ) R g , respectively. 
Assuming a 2-10 keV flux in 2002 of 3.9 x 10" 11 erg cur 2 
s _1 from RXTE-PCA monitoring, a luminosity distance 
of 41.3 Mpc (following Mould et al. 2000, and using H Q = 
70 km s- 1 Mpc- 1 and A Q = 0.73), the 2-10 keV luminos- 
ity £2-10 = 7x 10 42 erg s _1 . From Marconi ct al. (2004), 
an AGN with this £2-10 has a bolometric luminosity 
Lbo\ = 15^2-10 = 1-1 x 10 44 erg s _1 . The accretion rela- 
tive to Eddington m = £boi/£Edd is thus estimated to be 
0.15 for the 2002 XMM-Newton observation, m during 
the Suzaku observation is thus 0.04. Meanwhile, Lu & 
Wang (2000) have derived m ~ 5.5% from SED fitting. 
These values of m and our derived value of r t are not 
immediately consistent with the rough relation of Yuan 
& Narayan (2004), thus pointing toward models incor- 
porating smaller values of r t (see below). Furthermore, 
most low-m sources are radio loud, but NGC 4593 is not 
strongly radio-loud. Its 5 GHz flux has been measured to 
near ~2 mJy (e.g., Schmitt et al. 2001), and its B-band 
flux is ~6-16 mJy (e.g., McAlary et al. 1983), so the ra- 
dio loudness parameter, defined as the ratio of these two 
values, is <3. Values >10 define a source as radio-loud 
(Kellermann et al. 1989). A connection between rh (pro- 
portional to the observed X-ray flux) and r t in NGC 4593 
is thus qualitative only as well as speculative, especially 
since we have only two model-dependent estimates of r t . 

There is also the question of whether the inner portions 
of a thin disk in AGN can evaporate and/or become ra- 
diatively inefficient on timcscalcs of only a few years. As 
the accretion disks of BH XRBs are thought to evolve 
on timescales of at least hours to days, the correspond- 
ing timescales in NGC 4593 (black hole mass a factor of 
10~ 5 higher) would be decades to centuries. On the other 
hand, Marscher et al. (2002) interpreted rapid (duration 
of ^ a couple weeks) dips in the X-ray light curve of the 
radio-loud Seyfert 3C 120 as periods when the inner por- 
tion of the disk evaporated, each event leading to ejection 
of material into the relativistic jet and a corresponding 
radio flare about a month later. 

The total Fe Ka line intensity decreased from 4.79 ± 
0.48 x 10~ 5 (2002) to 2.93 ±0.23 x 10~ 5 (2007) ph cm~ 2 
s _1 , a factor of 1.7, while the observed 4-10 keV flux de- 
creased from 2.47 x lO" 11 (2002) to 0.69 x lO" 11 (2007) 
erg cm -2 s _1 , a factor of 3.6, i.e., roughly half of the total 
line flux has responded to continuum decrease. Modeling 
the Fe Ka line with a dual-Gaussian model attempted to 
separate the variable and non-variable emission compo- 
nents; in this model, a non- variable, narrow component 
is detected in both observations, and dominates the total 
line flux in the Suzaku spectrum, while a broader com- 
ponent is detected only in the XMM-Newton spectrum. 
The best-fit width a of the broad line was 

177^2 eV, 

corresponding to FWHM velocity of 19100±5goo km s _1 , 
implying r = 3.2±f f x 10 12 m = 0.12±%$. lt-days, or 

330^180-Rg- This estimate is inconsistent with B07's es- 
timate of ^2000 R g based on the invariance of the Fe 
Ka line during the XMM-Newton observation, but it 
is still consistent with the presence of a truncated thin 
disk (r t > 150i? g ). r t could be invariant from 2002 to 
2007, but an annular region on the thin disk spanning 



from ~ 300i? g to 1000-5000 R g (outer radius obviously 
speculative) could be responding to the drop in illumi- 
nating continuum flux. If the drop in Fe XXVI flux is 
real, then that line could also originate in this region. 
However, it is not clear in this model why material at 
> 6000i? g (yielding the narrow line component), con- 
tributing roughly half of the total line intensity in 2002, 
has not responded to the drop in continuum flux, as it is 
well within a week's light-travel time. The width of the 
narrow line had been fixed at a — 41 eV in our model- 
ing of the XMM-Newton profile, but contributions from 
a more narrow component likely cannot be ruled out. 
Such distant material could be responding to a previ- 
ous higher continuum flux. The 2-10 keV RXTE-PCA 
monitoring light curve in fact showed a higher, more av- 
erage flux level until roughly 300 days before the Suzaku 
observation (Figure 7). 

A final possibility that does not require evolution in r t 
is that the inner disk may have become be too ionized 
to transmit an Fe line. In the context of models with a 
hot, ionized skin (Nayakshin, Kazanas & Kallman 2000), 
a disk illuminated by a power-law continuum with a pho- 
ton index near 1.6, similar to that in the Suzaku obser- 
vation, yields extremely weak Fe line emission. 

Of course, both profile models are likely oversimplifica- 
tions. The community could thus benefit from an X-ray 
observatory with ~ few eV resolution combined with a 
large effective area near 6 keV to resolve the various com- 
ponents of the Fe Ka core as a function of time and/or 
continuum flux level, if multiple components do indeed 
exist, as well as resolve the Fe XXVI line. 

5.3. The Compton Reflection Component 

A Compton shoulder was not significantly detected 
in either the Suzaku or XMM-Newton spectra; we 
find upper limits to Compton shoulder emission (first- 
scattering) of ^23% of the core. It is thus not obvious 
from this limit alone whether bulk of the Fe Ka line orig- 
inates in Compton-thick material, especially the degree 
to which the strength of the Compton shoulder depends 
on the geometry of the material. However, no relativisti- 
cally broadened Fe Ka line has been confirmed in NGC 
4593, but we confirm from the broadband Suzaku spec- 
trum the presence of a Compton reflection hump which 
thus must correspond to (at least some fraction of) the 
Fe Ka core emission. We can investigate if the measured 
strength R of the Compton reflection hump, 1.08 ± 0.20 
(statistical uncertainty only; ± 0.35 including the sys- 
tematic uncertainty), can correspond to the observed Fe 
Ka line EW of 255±19 eV. Following George & Fabian 
(1991), one expects an EW of 135 eV (using the abun- 
dances of Lodders 2003) to correspond to R = 1 for an 
semi-infinite optically thick slab illuminated by a power- 
law continuum with T = 1.7 and assuming solar abun- 
dances and an inclination angle of 30° relative to the 
observer's line of sight. The observed values of R and 
EW are thus consistent if the Fe abundance relative to 
solar, Zp e , is about 1.7, which is not unreasonable. For 
a truncated disk, this could be explained by having the 
Comptonizing corona consist of numerous flares lying in 
a sandwich-like geometry just above/below the thin disk 
(e.g., Haardt et al. 1994), such that the disk spans 2ir sr 
of the sky as seen by each X-ray flare. However, very dis- 
tant (pc-scale), Compton-thick material lying out of the 
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line of sight, which cannot be ruled out as contributing to 
the observed Fe K emission profile, could also contribute 
to the total Compton reflection strength. 

However, if the thin disk is truncated, then a semi- 
infinite slab may not be an appropriate geometry, espe- 
cially if the central X-ray source is not immediately close 
to the reflecting disk. The EW of a truncated disk will 
be lower, but will depend on the location of the illumi- 
nating X-ray source. If we assume that the illuminat- 
ing X-ray source is located on the disk symmetry axis a 
height h above the disk, we can use the EW to constrain 
h. George & Fabian (1991, their Fig. 15) demonstrate 
that the reflected flux is dominated by the region of the 
disk with r/h ~ 1-2. For a truncated thin disk with r t = 
several thousand R gl h must also be ~ several thousand 
R g above the black hole. A configuration in which the 
X-ray corona is associated with the base of a jet along 
the symmetry axis may thus be applicable, e.g., Markoff, 
Nowak & Wilms (2005). NGC 4593, like many Seyferts, 
is known to host a pc-scale radio component (size of <15 
pc; Schmitt et al. 2001). 

5.4. Spectral Variability of the Broadband Components 

The primary-law component in Seyfert X-ray spectra 
is usually attributed to inverse Comptonization of soft 
seed photons. In an ADAF flow, thermal Comptoniza- 
tion is expected to dominate the X-ray emission unless 
the accretion rate is extremely low, in which case thermal 
bremsstrahlung emission dominates the X-ray spectrum 
(Narayan et al. 1998, Narayan 2005). One of our main 
results is that while the primary power-law component 
has dropped in flux by a factor of almost 4, the soft ex- 
cess has dropped in flux by a factor of ~20 between 2002 
and 2007, ruling out an origin for the soft excess with a 
size greater than 5 lt-years. This drop is likely linked to 



the decrease in the primary X-ray power-law, i.e., it may 
be either a cause of an effect of it. 

We explored two phenomenological models for the soft 
excess, bremsstrahlung and thermal Comptonization. In 
the presence of an ADAF flow, one can expect thermal 
bremsstrahlung emission with a temperature of 10 9 ~ n K 
(Narayan & Yi 1995), but such temperatures are higher 
by over 2 orders of magnitude compared to the tem- 
perature derived from our model fits and by B07. We 
also modeled the soft excess as inverse Comptonization 
of soft seed photons with an assumed input temperature 
of 50 eV by an optically thin corona with a tempera- 
ture k^T ~ 50 keV, again obtaining similar results to 
B07. The location of such a process is not clear, though 
it could occur in the ionized skin of the thin disk (e.g., 
Magdziarz et al. 1998, Janiuk et al. 2001), or at the base 
of an outflowing jet. If both the soft excess and hard 
X-ray power-law components originate via Comptoniza- 
tion of disk seed photons, a decrease in the intensity of 
those soft seed photons (e.g., from an increase in rt) be- 
tween 2002 and 2007 could yield a correlated drop in 
both component's flux. Another possibility is that the 
optical depth of the Comptonizing components may have 
changed. 
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APPENDIX 

In this Appendix, we discuss the X-ray spectrum of 1WGA J1239. 7-0526, an X-ray source located about 5' south 
of NGC 4593, at a = 12h39m42.8s, 5 = -5d26ml5s (J2000.0 coordinates from the NASA Extragalactic Database), 
and in the Suzaku-XIS field of view during the observation of NGC 4593. Detection of this source in X-rays was first 
reported by ROSAT-PSPC (White, Giommi & Angelini 2000). Its Galactic latitude is +57°, so it is likely not Galactic 
in nature. Its redshift is unknown. As this source is also in the field of view of the (non- imaging) HXD-PIN observation 
of NGC 4593, we extracted spectra of 1WGA J1239. 7-0526 from the XIS to estimate the level of contamination in the 
HXD-PIN spectrum of NGC 4593. 

Spectra were extracted in a manner similar to Section 2.1. ARFs were generated considering the source's position 
5' off-axis. Spectra were grouped to a minimum of 20 ct bin -1 . 

We fit the co-added 0.4-11 keV XIS-FI and 0.3-9 keV XIS-BI spectra in XSPEC with a model consisting of a simple 
power-law, with absorption by the Galactic column (1.88 x 10 20 cm~ 2 , Kalbcrla et al. 2005). This model yields a good 
fit with excellent data/model residuals, with \ 2 /dof = 369.0/378 for T = 1.75 ± 0.05 and a power-law normalization 
of 1.05 ±0.05 x 10~ 4 ph cm" 2 s" 1 keV" 1 at 1 keV. Absorbed fluxes for the 0.4-2 and 2-10 keV bands are 2.4 x 10" 13 
and 3.9 x 10~ 13 erg cm~ 2 s _1 , respectively. Background fluxes were 0.31 (FI) and 0.53 (BI) times the modeled source 
flux in the 0.4-2.0 keV band, and 0.58 (FI) and 2.2 (BI) the modeled source flux in the 2-10 keV band. 

Extrapolating this model to the HXD-PIN band yields a 12-76 keV flux estimate of 7.3 x 10~ 13 erg cm~ 2 s _1 , or 
about 1.5% of the net source flux of NGC 4593 in the 12-76 keV band. However, the HXD aimpoint is about 4' 
away from 1WGA J1239. 7-0526; at this angular distance, the PIN effective area is about 90% of that on-axis; 1WGA 
J1239.7-0526 thus contributes only ~1.35% of the net count rate of NGC 4593 in the PIN. 

To test for the presence of variability in 1WGA J1239. 7-0526, we extracted light curves from each XIS, binned them 
to 34560 s (6 satellite orbits), and added light curves from all three XlSes. In the 0.4-2.0 keV band, an upper limit to 
F val of 15.2% is found. In the 2-10 keV band, we measured the variability amplitude F var = 9 ± 6%. 
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TABLE 1 
Fe K Bandpass Model Parameters 







2007 Suzaku 


2002 XMM-Newton 


Component 


Parameter 


XIS-FI value 


EPIC-pn value 



Model Flux 
Power-law 
Compton Reflection 
Fc Ka emission line 



Fc XXV emission 



Fc XXVI emission line 



Compton Shoulder emission 



X 2 /dof 

4-10 keV flux (erg cm 2 

r 

R 

Energy (keV) 

Width a (eV) 

Intensity (ph cm -2 s _1 ) 

EW (eV) 

Energy (keV) 

Width a (eV) 

Intensity (ph cm~ 2 s _1 ) 

EW (eV) 

Energy (keV) 

Width a (eV) 

Intensity (ph cm~ 2 s _1 ) 

EW (eV) 

F-test significance 

Intensity (ph cm~ 2 s _1 ) 

EW (eV) 



320.6/318 

0.69 x 10- 11 

1.647 ± 0.043 

1.08 (fixed) 

6.412±0.007 
41 +12 

^-16 

2.93 ± 0.23 x 10" 5 
255±19 
6.70 (fixed) 
1 (fixed) 

< 2.8 x 10~ 6 

< 19 
6.96 (fixed) 

1 (fixed) 

< 1.8 X 10~ 6 

< 17 
N/A 

< 7 x 10" 6 

< 44 



1105.0/1091 
2.47 x lO" 11 

1.08 (fixed) 
6.396 ± 0.012 
87 ± 17 
4.79 ± 0.48 x 10" s 
114 ± 11 
6.70 (fixed) 
1 (fixed) 

< 5.5 x 10- 6 

< 12 
6.96 (fixed) 

1 (fixed) 
4.4 ±2.8 X 10~ 6 
12 ± 8 
92% (97.5%) 

< 1.1 x 10~ 5 

< 22 



Note. — Results are shown for the respective best-fit models for the 4-11.5 keV 2007 Suzaku XIS-FI spectrum and for the 4-11.5 keV 
2002 XMM-Newton pn spectrum in which the Fe Ka line is modeled as a single Gaussian. In each case, the best-fit model consists of a 
power-law component, a Compton reflection component (pexrav), and emission lines due to Fe I Ka and Fe I K/3; a Fe XXVI emission line 
was included in the fit to the XMM-Newton pn data. In the Suzaku data, the Fc Ka line energy has been corrected to account for the gain 
shift of +9 eV as measured using the 55 Fe calibration source spectra. The row labelled "F-test significance" denotes the confidence level 
associated with including this component in the model according to an F-test (omitted where line intensity is an upper limit only); values 
in parentheses denote statistical significance derived from Monte Carlo simulations (see Section 3 for details). 



TABLE 2 

Suzaku Broadband Spectral Fit Parameters 



Component 



Parameter 



Value 



Model Flux 
Power-law 

Warm Absorber 

Compton Reflection 
Brcmss. Component 



Fe L 3 Edge 



O VII / Emission Line 



Ne IX / Emission Line 



Mg XI / Emission Line 



XIS-BI/XIS-FI 



X 2 /dof 

2-10 keV flux (erg cm" 2 s" 1 ) 

r 

Norm, at 1 keV (ph cm~ 2 s _1 keV -1 ) 

Column density Nu.wA (cm -2 ) 

log(£) (erg cm s -1 ) 

Refl. strength R 

k B T (keV) 

Norm. 

Unabsorbed (absorbed) F .4_2.0keV 
Energy (keV) 



F-test significance 
Energy (keV) 
Intensity (ph cm~ 2 
EW (eV) 
F-test significance 
Energy (keV) 
Intensity (ph cm~ 2 
EW (eV) 
F-test significance 
Energy (keV) 
Intensity (ph cm~ 2 
EW (eV) 
F-test significance 
Norm. 



- 1 ) 



- 1 ) 



- 1 ) 



1254.5/1106 
1.03 x 10- 11 
1 656+ ' 018 
2.08 ±0.04 x 10~ 3 

2.8 ±0.8 x 10 21 
2.39 ± 0.17 
1.08 ± 0.20 

20+ 06 
0.8±q 5 x 10- 3 
^a'lt^xlO- 
0.707 (fixed) 
09+ 03 

u ' ua -0.04 

99.7% (96.4%) 
0.561 (fixed) 

2.9 ± 2.8 x 10" 6 
4.9 ±4.7 

90% (84%) 
0.905 (fixed) 

x lo- 5 



4.5" 1 



5 ±2 3 

98.4% (95%) 
1.331 (fixed) 



4.0 



+3.1 
1.6 

3.0 



x 10 

+2.3 



-6 



98.9% (94%) 
1 059+ ' 010 



Note. — Best-fit model parameters for the Suzaku broadband fit with the soft excess modeled as thermal brcmsstrahlung emission. See 
the XSPEC user manual for the units of the bremsstrahlung normalization. The rows labeled "F-test significance" denote the confidence 
level associated with including this component in the model according to an F-test; values in parentheses denote statistical significance 
derived from Monte Carlo simulations (sec Section 3 for details). The uncertainty on R is statistical only; the systematic uncertainty is an 
additional ±0.15. 
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TABLE 3 

Comparison of Broadband Model Fits to XMM-Newton and Suzaku spectra 



Component 


Parameter 


2002 XMM-Newton 
EPIC-pn value 


2007 Suzaku 
value 


Thermal Bremsstrahlung model 


(Absorbed) Model Flux 

Power-law 

Warm Absorber 
Fe L3 edge 
Brcmss. Component 


2-10 keV flux (erg cm~ 2 s" 1 ) 
0.5-2.0 keV flux (erg cm" 2 s" 1 ) 

r 

Norm, at 1 keV (ph cm" 2 s" 1 keV -1 ) 
Column density -/Vjj,wa (cm -2 ) 
Optical depth tf c l3 
fe B T (keV) 
Norm. 

Unabsorbed (absorbed) -F0.4-2.okcV 


1793.23/1756 
3.93 x 10- 11 
2.57 x 10- 11 
1.901±0.006 
1.19 ±0.01 x 10- 3 
1.7 ±0.2 x 10 21 
0.16 ±0.2 
0.28 ± 0.01 
2.02 ±0.01 x 10~ 2 
4.72 ±0.19 (3.87 ±0.15) xlO" 12 


1254.8/1105 

1.03 x 10 -11 

4.15 x 10~ 12 

1 656+ - 018 
i.dod_ 025 

2.08 ±0.04 X 10- 3 

2.8 ±0.8 x 10 21 

o.o9tg:S 
o.2ot8:« 

0.8+Q5 x 10_3 
S.VtlY (7.lt 5 ; 9 ) X10-" 


Thermal Comptonization model 



(Absorbed) Model Flux 
Power-law 

Warm Absorber 
Fe L3 edge 

Comptonized Component 



2.0 keV flux (erg cm" 2 s" 1 ) 



7, 
0.5 
r 

Norm, at 1 keV (ph cm" 2 s~ 
Column density -/Vjj.wa (cm' 
Optical depth rp c L3 
Temperature (eV) 
Optical depth r 
Norm. 

Unabsorbed (absorbed) F0.4. 



1 keV 

- 2 ) 



- 1 ) 



1945.8/1757 
2.57 x lO" 11 
1 75n +0 ' 018 

8.30+.?;? 9 , x 10_: 
2 i 3+ - 3 ^ in2l 



1255.775 / 1105 
4.15 x 10~ 12 



-0.2 
0.15 4 



x 10 2 
ho.01 
-0.02 
53 (fixed) 
n n +0.03 
U -' i -0.05 

i-GOtoH x 10" 



2.09 
3.1 



1 65Q+ - 13 
l.DO»_ _ 011 

+0.05 

0.04 



x 10" 



-1.2 
-0.9 



x 10 



21 



1.54 



+0.33 
0.13 



(l.28l8;£) x!0" 



8.2 



0.08±0.04 

53+ 57 
<0.28 
2.8t 2 8 3 7 x 10~ 5 
+54 7 6 (6.7l 5 4 5 4 ) xlO" 



Note. — Results from applying the best-fit broadband models from the Suzaku observation to the 0.2—12 keV XMM-Newton pn 
spectrum, using thermal bremsstrahlung or thermal Comptonization components to model the soft excess. Col. (2) lists the parameters 
which were allowed to vary from their values in the Suzaku fit; the values for the XMM-Newton fits are listed in Col. (3). The corresponding 
values from the Suzaku fits are listed in Col. (4). See Table 2 for other model parameters for the thermal bremsstrahlung fit. Model fluxes 
are in units of erg cm" 2 s _1 . Sec the XSPEC user manual for the units of the Comptonized component. 




Fig. 1. — Binned light curves from the Suzaku observation. The top panel shows the 0.4-10 keV light curve, summed over XIS 0+1+3 
and binned to 5760 s. The bottom panels shows the 12-76 keV PIN light curve, binned to 11520 s; error bars include a ~2% systematic 
uncertainty (Fukazawa ct al. 2009). The x-axis denotes time since 2007 December 15 at 04:09 UT. 
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Fig. 2. — Data/model residuals to fits to the Suzaku XIS-FI Fc K bandpass spectrum. Data have been rebinned by a factor of 4. Residuals 
to a simple power-law fit are shown in panel (a). Panel (b) shows residuals when the Fe Ka line is modeled with a Gaussian. In panel (c), 
a Gaussian to model Fe Kf3 line emission has been included. In panel (d), a Compton reflection component, containing an Fe K edge at 
7.11 keV, has been included; this is our best-fit model. Vertical dashed lines denote the expected energies of Fe I Ka, Fe XXV, Fe XXVI, 
and Fc I K/3 emission lines. 
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Fig. 3. — Data/model residuals to fits to the Suzaku XIS-FI (black points; rebinned by a factor of 7) and 2002 XMM-Newton pn (gray 
points; rebinned by a factor of 9) Fe K bandpass spectra. Vertical dashed lines denote the expected energies of Fe K a, Fe XXV, Fe XXVI, 
and Fe K/3 emission lines. Panel (a) shows residuals to a model consisting of a power-law plus a Compton reflection hump; panel (b) shows 
residuals when the Fe Ka and K/3 lines are included in the model as well. 



I 

£ 

o 
ft 



1 

0.5 


-0.5 



1 1 1 

- 2007 Suzaku 

: 


- 










^-f— 1 1 


— 1 1 1 — F 


- 2002 XMM pn 


- 




_ 



I 

£ 

o 

ft 




6.5 6.6 6.7 6.8 
Energy (keV) 



6.9 



Fig. 4. — Contour plots resulting from applying a "sliding Gaussian" to a model consisting of a power-law, a Compton reflection component, 
and the Fe Ka line. The resulting plots of line intensity versus energy are shown for the Suzaku XIS-FI (upper) and XMM-Newton pn 
(lower) Fe K bandpass spectra. The areas outside the contours indicate regions of (energy, intensity) parameter space where adding a 
Caussian to this model would make the fit worse The rest-frame energies of Fe XXV and Fe XXVI are denoted by vertical dashed lines. 
Solid, dashed, dot-dashed, and dotted lines denote Icr, 2a, 3cr, and 4<r confidence levels, respectively, for two interesting parameters. Note 
that the y-axis scales are the same. This figure indicates that the Fe XXVI line detected in the XMM-Newton spectrum is not significantly 
detected in the Suzaku spectrum. 
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Fig. 5. — Data/model residuals to the broadband spectrum for the Suzaku XIS-FI (black points < 12 keV), XIS-BI (grey points), and 
HXD-PIN (black points > 12 keV). The data are plotted in panel (a). Residuals to a simple power-law fit (panel (b); assuming PIN/X1S-F1 
normalization of 1.16) show the moderately strong soft excess as well as Fe L absorption features near 0.9-1.2 keV associated with the 
warm absorber. The hard excess above 10 keV due to the Compton reflection hump is also evident. In panel (c), the Compton reflection 
component, warm absorber, Fe K emission lines and soft excess (modeled as thermal brcmsstrahlung emission) have been modeled. In 
panel (d), the Fe L3 edge and the He-like lines of O, Ne, and Mg have been added. Vertical dashed lines denote the expected energies of 
the O VII line, the Fe L3 edge energy, the Ne IX line, and the Mg XI line. All data are rebinned by a factor of 4. 
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Fig. 6. — Contour plot of Compton reflection component strength R as a function of power-law photon index T, assuming the PIN/XIS-FI 
normalization is fixed at 1.16. The cross marks the best-fit parameters. Solid, dashed, dot-dashed, and dotted lines denote lc, 2<r, 3cr and 
4cr confidence levels, respectively, for two interesting parameters. 
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dashed line denotes the long-term average mean flux of 3.9 X 10 
XMM-Newton and 2007 Suzaku observations 
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Fig. 8. — Data/model residuals to our best-fit to the 2002 XMM-Newton pn spectrum. Data are rebinned by a factor of 8. 
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Fig. 9. — Unfolded model spectrum illustrating the best-fit model to the Suzaku spectrum. Red crosses denote the data, rebinned by a 
factor of 4. Plotted are the power-law component (red), Compton reflection component (green), thermal bremsstrahlung component (blue), 
He-like O, Ne and Mg lines (yellow), and Fe Ka and K/3 emission lines (cyan). 
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Fig. 10. — Unfolded model spectrum illustrating our best-fit model to the 2002 XMM-Newton pn spectrum. Data are binned by a factor 
of 6. Components are denoted by the same colors as in Fig. 9; the Fe XXVI emission line is also plotted in orange. 



